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20.  Abstract  (Continued) 

)  previous  study  of  satellite  based  radiances.  Possible  sources  of  this  dis¬ 
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An  Atmospheric  Temperature  Profile 
Measured  With  an  In-Situ  Infrared  Radiometer 


1.  INTRODl  CnON 

A  great  deal  of  effort  has  been  expended  during  the  past  15  years  to  develop 
hardware  and  analysis  techniques  to  derive  vertical  distributions  of  atmospheric 
temperatures  from  radiance  measurements  made  from  a  satellite-based  sensor. 
Most  of  these  measurements  have  made  use  of  the  15  *jm  carbon  dioxide  band  as 
first  suggested  by  Kaplan^  and  have  been  based  on  the  fact  that  carbon  dioxide  can 
be  taken  as  uniformly  mixed  in  the  free  atmosphere.  This  assumption  allows  us 
to  make  a  rather  good  estimate  of  the  transmission  properties  of  the  atmosphere 
in  the  spectral  region  of  interest.  In  principle,  then,  the  radiances,  as  meas¬ 
ured  in  different  spectral  intervals  of  the  15  ym  carbon  dioxide  band,  ranging 
from  the  band  center  to  one  of  the  band  wings,  can  be  interpreted  in  terms  of  the 
vertical  distribution  of  atmospheric  temperature.  Figure  1  indicates  the  basic 
concept  calculated  for  the  U.S.  Standard  Atmosphere,  1976,  in  terms  of  the 
spectral  channels  that  compose  the  temperature  sounder  on  the  Defense  Meteoro¬ 
logical  Satellite. 


(Received  for  publication  1  May  1979) 

1.  Kaplan,  L.  D.  (1959)  Inference  of  atmospheric  structure  from  remote  radiation 

measurements,  J.  Op.  Soc.  Amer.  49:1004. 

2.  U.S.  St^dard  Atmosphere  (1976)  NOAA-S/T  76-1562.  Supt.  of  Documents, 

U.  S.  Government  Printing  Office. 


Fijfure  1.  Calculated  Spectral  Distribution  of  I'pwelling  Radiance  in  the 
15  pm  Carbon  Dioxide  Band.  Indicated  pressure  levels  show  the  portion  of 
the  atmosphere  responsible  for  the  maximum  upwelling  emission  to  space 
at  the  frequencies  of  the  DMSP  filters 


Since  the  concept  of  remote  temperature  sounding  of  the  atmosphere  was  first 
suggested,  the  development  of  objective  inversion  algorithms  has  been  fraught 
with  difficulty  and  largely  discarded  in  favor  of  statistical  techniques  requiring 
climatological  data  as  a  basis  of  comparison.  In  an  effort  to  get  at  the  root  of 
this  difficulty,  a  project  was  designed  to  investigate  and  verify  our  ability  to  solve 
the  "forward  problem"  that  is,  the  calculation  of  the  expected  spectral  radiances 
emerging  from  an  atmosphere  whose  profile  of  temperature  and  moisture  is  known 
from  radlosonde/rocketsonde  measurements. 

The  Defense  Meteorological  Satellite  Program  (DMSP)  includes  a  sensor  cap¬ 
able  of  measuring  radiances  in  the  15  pm  region  at  frequencies  located  as  indicated 
in  Figure  1.  Using  several  sets  of  measurements  from  the  DMSP  satellite,  a 
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a  atudy  of  the  "forward  problem"  was  conducted  by  McClatchey,  Thla  study 
began  with  the  identification  of  a  number  of  satellite  measurements  made  In  geo¬ 
graphical  regions  in  close  proximity  in  time  and  space  to  rawinsonde/rocketsonde 
stations,  h'urther,  an  effort  was  made  by  inspection  of  visual  and  Infrared  Imagery 
obtained  by  sensors  on  the  same  satellite,  to  establish  that  the  atmospheric  path 
from  surface  to  satellite  was  free  of  clouds.  Carbon  dioxide  was  assumed  uni¬ 
formly  mixed  and  the  rawinsonde/rocketsonde  temperature  data  were  used  to  com¬ 
pute  the  upwelling  radiation  corresponding  to  the  six  filter  channels  of  the  satel¬ 
lite  sensor.  With  the  exception  of  one  channel  located  near  the  carbon  dioxide 
band  center,  all  channels  for  all  cases  studied  showed  a  significam  lysfcrnatic 
discrepancy,  the  calculated  values  being  larger  than  those  measured.  The  mag¬ 
nitude  of  this  discrepancy  is  largest  in  the  707  cm*^  channel  (15  percent)  and 
smallest  in  the  Q-branch  channel  centered  at  668  cm  *. 

Three  potential  sources  of  these  radiance  discrepancies  have  been  considered; 
1)  A  sensor  calibration  error;  2)  An  uncertainty  in  the  calculated  transmittances; 
and  3)  Rawinsonde/rocketsonde  temperature  errors  (or  some  difference  between 
the  rawinsonde  tempei-atures  and  a  radiatively  derived  temperature).  The  ques¬ 
tion  of  instrument  calibration  has  been  investigated  and  potential  uncertainties 
due  to  this  cause  are  felt  to  be  far  less  than  those  observed.  The  transmission 
uncertainty  issue  is  being  investigated,  and  current  estimates  are  that  it,  too,  is 
unlikely  to  be  the  cause  of  the  large  radiance  discrepancy  noted.  The  third  item 
is  the  subject  of  this  report. 

2.  RARIOMETRiC  AND  .St;i>rORTING  METEOROLOGICAL 

MEASUREMENTS 

A  liquid-nitrogen-cooled  radiometer  operating  in  the  15  pm  region  was  flown 
on  a  balloon-borne  platform  launched  from  White  Sands,  New  Mexico  on 
28  September  1977  (see  Murcray,  et  al^).  The  Intent  of  this  measurement  was  to 
filter  the  radiometer  in  the  strongest  absorbing  portion  of  the  15  pm  carbon  diox¬ 
ide  band  and  to  view  the  atmosphere  horizontally.  Transmission  calculations 
indicated  that  the  atmosphere  would  be  opaque  over  short  horizontal  distances  all 
the  way  from  the  surface  to  an  altitude  of  20  km  and  that  even  at  30  km  (the  max¬ 
imum  height  of  the  balloon)  the  atmosphere  would  become  opaque  over  a  horizontal 
distance  of  a  few  hundred  kilometers.  Therefore,  In  the  lower  portion  of  the 
atmosphere,  the  radiometer  would  be  expected  to  measure  the  local  temperature. 

3.  McClatchey,  Robert  A.  (1976)  Satellite  Temperature  Sounding  of  the 

Atmosnhere;  Ground  Truth  Analysts,  AFuL-TR-76-0279.  AD  A038236. 


A.  higher  altitudes,  a  comparison  between  measurements  and  calculations  could 
be  made  via  the  radiative  transfer  equation  and  the  presumed  known  temperature 
profiles  measured  by  rawinsonde  and  rocketsonde. 


The  liquid  nitrogen  used  to  cool  the  entrance  aperture  of  the  radiometer  is 
allowed  to  vent  out  through  an  entrance  baffle  system  which  both  cools  the  baffle 
and  provides  an  antifrost  capability  for  the  entrance  window.  In  order  to  avoid 
dumping  liquid  nitrogen,  it  was  necessary  to  orient  the  radiometer  at  an  elevation 
angle  of  20  .  This  reduces  somewhat  the  altitude  at  which  the  radiometer  can  be 
considered  to  be  monitoring  the  local  temperature.  However,  this  situation  can 
countered  bv  applying  a  radiative  transfer  equation  at  somewhat  lower  altitudes 
than  originalh  intended.  A  high  resolution  spectrum  based  on  the  f.  S.  Standard 
Atmosphere.  1976  and  the  actual  measurement  geometry  is  given  in  Figure  2. 

This  figure  helps  provide  an  understanding  of  the  actual  atmosphere  opacity  at 
float  altitude. 


The  radiometer  was  carefully  calibrated  against  a  known  black  body  on  four 
separate  occasions.  No  difficulties  were  encountered  and  the  four  calibration 
results  were  consistent  to  better  than  0.  5  percent.  Consideration  of  other  possi¬ 
ble  errors  associated  with  calibration  procedures  leads  to  the  confident  statement 
that  the  radiance  measurements  should  not  be  as  large  as  1  percent 

The  radiometer  was  filtered  so  that  the  net  spectral  response  function  of  the 
instrument  is  as  indicated  in  Figure  3.  Although  it  was  originally  intended  to  use 
a  narrow  filter  located  in  the  Q-branch  (see  Figure  2)  (about  1.  5  cm'*  wide  cen¬ 
tered  near  668  cm  ).  it  was  not  possible  to  obtain  such  a  filter  without  great 
expense.  Therefore,  an  available  filter  centered  near  683  cm'^  near  the  strong¬ 
est  portion  of  the  R -branch  was  used  instead.  This  fact  also  contributes  to  a 
decreased  opacity  at  higher  altitudes.  This  lower  opacity  simply  requires  that 
the  measurements  be  analyzed  by  application  of  the  radiative  transfer  equation  at 
somewhat  lower  altitude  than  would  otherwise  be  required.  Figure  4  has  been 
included  to  provide  an  indication  of  the  net  atmospheric  opacity  as  viewed  from 
the  radiometer.  It  can  be  seen  that  the  atmosphere  becomes  significantly  non¬ 
opaque  above  150  mb. 


Figure  5  shows  the  details  of  the  flight  including  launch  time,  the  time  the 
sensor  was  located  at  certain  pressure  altitudes,  the  general  path  of  the  balloon 
package,  and  the  location  of  rawinsonde  ground  stations  from  which  vertical  pro¬ 
file  information  was  obtained.  Figure  6  shows  the  temperature  field  at  the  30  mb 
pressure  level  as  obtained  from  rawinsonde  information.  This  figure  is  shown  as 
an  indication  of  the  temperature  uniformity  over  a  broad  geographical  region.  At 
this  pressure  altitude,  the  temperature  can  be  seen  to  be  uniform  to  within  1  or 
2  C  over  distances  of  3000  km.  Temperature  data  at  lower  altitudes  show  a 
Similar  degree  of  uniformity.  This  fact  is  important  as  we  use  the  local  rawinsonde 
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Figure  2.  Calculated  High  Resolution  Spectrum  of  Car¬ 
bon  Dioxide  for  a  U.  S.  Standard  Atmosphere,  1976  (Ref. 
2),  Looking  Upward  at  a  20°  Elevation  Angle  at  an  Alti¬ 
tude  of  30  km 


FILTER  TRANSniSSION 


Figure  3.  The  Transmission  Function  of  the  Filter  Used 
in  the  Side-Looking  Radiometer 


Figure  4,  The  Transmission  From  the  Indicated  Pressure  Level  to  Space  ^ilong  a  20°  Elevation 
Angle  (70°  Zenith  Angle) 


Figure  5.  The  Ground  Track  of  the  Balloon  Payload  at  White  Sands  Missile 
Range,  New  Mexico 
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Figure  6.  The  30  mb  Temperature  Field  Over  Much  of  the  Northern 
Hemisphere  on  28  September  1978 


L  ATI  TUDE 


data  as  Input  to  the  radiative  transfer  equation  under  conditions  when  the  atmos¬ 
phere  becomes  more  transparent  at  the  higher  altitudes. 

^  igure  7  is  a  typical  example  of  one  of  the  rawinsonde  and  rocketsonde  tem¬ 
perature  profiles  above  the  200  mb  level.  The  rawinsonde  data  shown  in  this 
figure  correspond  to  the  0900  MDT  rawinsonde  launched  on  28  September  1977. 
The  rocketsonde  data  above  25  mb  can  be  seen  to  differ  by  a  few  degrees  from 
the  rawinsonde  data  in  the  region  of  overlap.  The  rocketsonde  was  launched  at 
1427  MDT  on  28  September  1977.  In  the  region  of  overlap  of  rawinsonde  and 
rocketsonde  data,  some  teats  were  run  to  study  the  effect  of  using  first  one  and 
*h;‘n  the  other  set  of  data  on  the  computed  upwelling  radiances.  Figure  8  shows 
the  result  of  tins  analysis  for  the  0900  rawinsonde  data.  In  view  of  the  small 
change  and  the  direction  o'  that  change,*  it  was  decided  to  use  the  rawinsonde  data 
to  the  maximum  possible  height  and  to  add  the  rocketsonde  data  at  the  top  of  the 
rawinsonde  sounding.  The  atmospheric  profiles  of  temperature  and  composition 
derived  from  rawinsonde  and  rocketsonde  data  used  in  all  calculations  described 
in  the  next  section  are  provided  in  Tables  la  to  If. 

3.  ANAI.Y.SI.S  OK  KAUIANCK  MEAStRKMKNTS 

Due  to  a  recovery  failure,  the  lower  altitude  data  (for  pressures  greater 
than  290  mb)  were  not  obtained.  However,  good  data  were  obtained  at  higher 
altitudes.  Based  on  the  atmospheric  opacity  as  indicated  in  Figure  4,  data  ob¬ 
tained  at  pressure  levels  greater  than  about  150  mb  might  be  expected  to  be 
close  to  the  local  atmospheric  temperature.  We  might  expect  the  deviation  in  the 
region  from  290  to  150  mb  to  increase  with  decreasing  pressure.  At  higher 
altitudes  (lower  pressures)  it  is  clear  that  the  radiances  received  at  the  sensor 
must  be  an  appropriate  weighted  average  of  radiances  along  the  atmospheric  path 
determined  by  application  of  the  radiative  transfer  equation.  Equation  (1)  is  the 
usual  form  of  the  transfer  equation  applicable  to  the  problem  assuming  the  atmos¬ 
phere  to  be  in  Local  Thermodynamic  Equilibrium  with  the  source  function  set 
equal  to  the  Planck  function. 

ijr  =  t)  dTo;)  (1) 


4i 

(This  will  tend  to  minimize  the 
In  Section  3. ) 


measurement -calculation  discrepancy  dt  cribed 
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Figure  7.  Typical  Rawlnaunde/Rocketaonde  Data 
Obtained  at  0900  MDT  on  28  September  1977 
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Figure  8.  Radiance  Profile  Computed  Based  on 
Rawlnaonde  Data  <^>  and  Rocketaonde  Data  (x)  In 
the  Region  of  Data  Overlap  as  Shown  In  Figure  7 
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Table  la.  Atmospheric  Profiles  of  Temperature  and  Composition 
Derived  from  Kawlnsonde  and  Rocketsonde  Data 


White  Sands  Missile  Range.  2tl  Sept. 

1977,  0900  MDT 
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39f779 

6.990C*28 

6.198E«17 

1.671E*ia 

ei.iM 

213.99 

6.290E629 

3.966E»17 

1.017E*10 

96. 999 

319.29 

9.961E«29 

3.966E»17 

1.979E*ia 

99.999 

319.65 

3.900Ee20 

2.979E»17 

I.297E*ia 

69.090 

319.29 

2.e90Ea20 

2.976E*17 

2.9992*18 

31.999 

319.19 

2.236E»29 

1.889E«17 

3.911E*10 

39.999 

321.65 

2.109E«20 

1.7a2E«17 

3.669E*10 

99.99f 

339.99 

?.916E«20 

1.711C« 17 

3.791E*18 

99.099 

331.69 

1.790Ee20 

1.699E»17 

«. 1702*10 

90.999 

329.99 

1.609Ee29 

1.19aEei7 

6.7072*10 

19.399 

399.79 

1.397E620 

1.191E«17 

6.869E*ia 

19.739 

331.99 

1.172E«'20 

9.963Eal9 

9.222E*18 

16.993 

332.05 

l.ai2Ea29 

9.591Eel6 

9.552E*18 

19.911 

336.99 

9.799E«19 

7.632Eei6 

9.a91E*18 

19.999 

396.99 

7.987E«1S 

6.63aEei9 

6.116E*ia 

9.1%7 

337.99 

9.7S3Eel9 

6.939E«16 

6.996E*ia 

9. 1ST 

362.35 

6.S10Eei9 

3.697Eei9 

C.991E*ia 

6.9S1 

359.95 

2.a96E»19 

2.639E*16 

3.3992*10 

9.197 

396.99 

1.696Eei9 

1.269E«16 

1.129E*17 

1.999 

379.99 

7.962Eai9 

5.976Eei5 

3.0712*16 

.939 

293.65 

3.796Eal9 

9.l99Eel5 

1.679E*16 

•  999 

35I.99 

2.296Eei9 

1.913E»19 

0.209E»19 

.993 

SS9.1S 

1.991E»16 

1.681E*19 

7.011E*19 

•  919 

399.69 

l.*19Eei8 

1.2a9E«19 

6.9992*15 

.199 

267.19 

1.330Eel8 

1.129Eel9 

6.175E*19 

Table  lb.  Atmoaphei-lc  Profiles  of  Temperature  and  Composition 
Derived  from  Rawlnsonde  and  Rocketsonde  Data 


White  Sands  Missile  Range,  2&  Sept. 

1977.  1250  MDT 

PRESSURE 

TERPERATURE 

MOLECULES/SC 

CP 

INt) 

IKl 

C02 

N20 

03 

SSS.Ett 

250.85 

2.519E»21 

2.137Eel8 

5.a67EP17 

351. lOf 

209.05 

2.050E421 

2.079Eel8 

5.130EP17 

3M.0IR 

201.35 

2.100Ee21 

1.702Eel0 

5.767E*17 

2«3.ai0 

237.75 

1.5a7E*21 

1.686E»18 

5.982C*17 

273. aot 

236.85 

1.917E»21 

1.626E*18 

6.122EP17 

ESt.OOO 

232.05 

l.T50E«21 

l.O05E«18 

€.053E*17 

210.000 

221.05 

l.OOOEVEl 

i.iaoEeia 

i.200E4l7 

17S.I00 

215.65 

1.225E»21 

l.OOOEeia 

7.633Eei7 

17k. 300 

215.05 

1.220E*21 

1.035Eei0 

7.606EP1 7 

ISO. 000 

209.95 

i.esoEeai 

8.910E*17 

8.180E«17 

lES.flOO 

202.15 

8.750E420 

7.025E*17 

I.958EP17 

12k.600 

202.05 

a.722£«20 

7.0aiE417 

6.973EP17 

ICO.OOO 

195.85 

7.C00E»20 

5.900E»17 

1.033EP10 

ao.ooo 

201.75 

5.600E«20 

0.752E»17 

1.251EeiO 

7a. 000 

205.35 

O.500E420 

O.150E«17 

1.071EP10 

eo.ooo 

211.55 

O.200E«20 

3.560Eel7 

1.017EP18 

56.300 

210.05 

3.901C»20 

3.300Eel7 

1.975E*16 

50.000 

210.65 

3.500E*20 

2.970E*17 

2.207EP18 

kO.OOO 

217.85 

2.600E>20 

2.376£«17 

2.890Eei0 

30.000 

222.05 

2.100E«20 

1.782E»17 

3.668EP18 

25.000 

223.85 

1.790E»20 

1.085E»17 

0.178E*18 

20.000 

226.05 

1.000E»20 

1.180E«17 

0.707EP18 

IS.SOO 

229.95 

1.106E^20 

9.385E416 

5.353E»iO 

15.000 

230.35 

l.e50C«20 

0.910E«16 

5.070E*18 

10.000 

233.65 

7.e00Eei9 

S.900E»16 

(.206EP18 

0.300 

233.05 

6.510E»19 

5.520Eei6 

f.359EP10 

a. 107 

237.95 

5.703E^19 

0.839E*16 

€.550E*18 

6.157 

202.35 

O.310E«19 

3.657Eei6 

C.901E«16 

k.OOl 

250.05 

2.660E«19 

2.O30E416 

3.35aEaia 

2.137 

260.65 

1.096Eei9 

1.269E»16 

1.128E*17 

1.006 

275.35 

7.CO2E«10 

5.976e»15 

3.071EP1E 

•  530 

263.05 

3.T66C«10 

3.196E«15 

1.079EP16 

.322 

258.05 

2.'250E«18 

1.913Eel5 

6.200EP15 

.203 

256.15 

i.saiE^io 

1.661E41S 

7.011EP15 

.217 

250.05 

1.5l9Cai0 

1.209Eei5 

0.999EP15 

.190 

207.15 

l.S30E*10 

1.129E«i5 

0.175E»15 

TftbU  tc.  Atmoapterlc  ProfU**  of  Tomporatur*  and  Compoattlon 
Darivad  from  Rawlnaonda  and  Rockataoada  Data 


Whlta  Sanda  Mtaalla  Ranfa,  28  Sapt 

.  1977,  1300  UOT 

rars5utc 

TfRPCRATURf 

MOiCctiLts/sa 

CN 

l««l 

(R) 

C02 

N20 

03 

!I8.888 

293.29 

2.C88Ca21 

2.2782*18 

8.7222*17 

282. 881 

291.89 

2.987Ca21 

2.1812*18 

8.9972*17 

198.188 

289.89 

2,898ta21 

2.8792*18 

9.1382*17 

288.888 

281,89 

2.188C«21 

1.7822*18 

9.7872*17 

288.888 

828.99 

1.882Ca21 

1.9972*18 

2.1932*17 

298.888 

222.99 

1.798Ca21 

1.8892*18 

9.8932*17 

288.888 

228.79 

t.88Ua21 

1.1882*18 

7.2882*17 

179.888 

218.89 

lJ229Ca21 

1.8882*18 

7.8132*17 

172.888 

212.79 

l.tl7Ca21 

1.8322*18 

7.8982*17 

198.888 

288.89 

t.8982a21 

8.9182*17 

8.1882*17 

129.888 

281.89 

8, 7982*28 

7.8292*17 

8.9982*17 

128.888 

881.f9 

8,7282*28 

7.8132*17 

8.9892*17 

189.888 

198.29 

7.8882*28 

8.2892*17 

9.9272*17 

188.888 

198.99 

7.8882*28 

9.9882*17 

1.8332*18 

91.888 

197.89 

8.2782*28 

9.8892*17 

1.1182*18 

81.888 

281.79 

9.8882*29 

8.7922*17 

1.2912*18 

78.888 

888.29 

8.98824»8 

8.1982*17 

1.8712*18 

88.888 

211.99 

8.2892*28 

3.9882*17 

1.8172*18 

98.888 

218.89 

3.8282*28 

3.2992*17 

2.8822*18 

98.888 

219.19 

3.9882*28 

2.9782*17 

8.2872*18 

88.888 

218.89 

2.8882*<28 

2.3782*17 

2.8982*18 

18.888 

221.79 

2.1882*28 

1.7822*17 

3.8802*18 

29.288 

229.89 

1.7712*28 

1.9832*17 

8.1882*18 

29.888 

229.79 

1.7982*28 

1.8892*17 

4.1782*18 

28.888 

227.89 

1.8882*28 

1.1882*17 

8.7872*18 

19.888 

229.89 

1.1982*28 

8.9182*18 

9.8782*18 

11.888 

'222.19 

8.2882*19 

7.8892*18 

9.9832*18 

18.828 

228.99 

7.9872*19 

8.8382*18 

9.1182*18 

8.1%7 

827,99 

9.7832*19 

8.8392*18 

8.9982*18 

8.197 

282.39 

8.818£*19 

3.8972*18 

9.9912*18 

8.891 

298.89 

2.8882*19 

2.8382*18 

3.2982*18 

2.127 

288.89 

1.8982*19 

1.2892*18 

1.1282*17 

1.888 

279.29 

7.8822*18 

9.9782*19 

3.8712*19 

.928 

183.89 

3.7882*18 

3.1982*19 

1.8792*18 

.322 

298.89 

2.2982*18 

1.9132*19 

8.2992*19 

.282 

298.19 

1.9812*18 

1.8112*19 

7.8112*19 

.217 

298.89 

1.9192*18 

1.2892*19 

8.9992*19 

.198 

287.19 

1.3382*18 

1.1292*19 

4.1792*19 

Tabl*  Id.  Atmoapharlc  ProfUea  of  Temper aturo  end  Compoaltlon 
Derived  from  Rewlnsoade  end  Rocketeonde  Dele 


White  Sends  Missile  Renge,  28  Sept.  1977.  2100  MDT 


RRfSStWC 

TtUUITtMC 

MClfCULlS/Se 

ce 

rail 

W1 

C02 

N20 

03 

Sll.lll 

169.69 

2.617E»11 

1.1666618 

6.9926 

S99.MI 

169.69 

8.691C6I1 

1.1136618 

9.8976 

Sll.tll 

139.89 

2.  luce'll 

1.7816618 

9.7676 

tIf.lM 

138.19 

l.U6Ce8l 

1.736C618 

9.8U6 

2IX.1H 

137.19 

1.976eefl 

1.6776618 

(.8826 

XSI.MI 

111.89 

1.798Ce21 

1.6996618 

(.6936 

Xll.lll 

117.69 

l.ClXCell 

1.377C618 

(.7126 

Xff.tll 

R1I.79 

l.U8C«*lt 

1.1886618 

7.2866 

tIS.SII 

119.99 

l.t83Ce21 

1.8896618 

7.6826 

ISI.MI 

187.69 

1.898Cell 

8.9186617 

(.1886 

tl6.su 

186.19 

9.9618688 

8.8966617 

(.9666 

itI.UI 

181.69 

8.6116628 

7.1666617 

9.1926 

117.9M 

186.19 

8.211C618 

6.9Ue6l7 

9.2916 

tlf.lll 

888.19 

7.6886618 

9.9686617 

1.8336 

U.3M 

197.79 

6.C71C618 

9.6616617 

1.8716 

I9.3M 

199.39 

9.8316618 

6.9686617 

1.2816 

71. Ill 

186.99 

6.9886618 

6.1986617 

1.6716 

M.UI 

119.99 

3.  9186618 

1.978C617 

2.2876 

II.SM 

116.19 

2.6816618 

1.1796617 

3.8176 

II.IM 

118.89 

1.1886618 

1.7816617 

3.(686 

II.HI 

116.19 

1.6886618 

1.1886617 

6.7876 

16.su 

126.89 

1.9816618 

8.6966616 

9.9776 

12.SU 

117.39 

8.7986619 

7.6296616 

9.8926 

11.131 

136.99 

7.9876619 

6.6386616 

(.1166 

1.167 

137.99 

9.7836619 

6.8396616 

(.9966 

6.117 

162.39 

6.3196619 

3.6976616 

(.9816 

6.191 

198.89 

1.8666619 

2.6386616 

3.3986 

1.137 

166.69 

1.6966619 

1.1(96616 

1.1X86 

I.IU 

179.39 

7.8616618 

9.9766619 

3.8716 

.9U 

163.69 

3.7666618 

3.1966619 

1.6796 

.321 

191.89 

2.2966618 

1.9136619 

8.2186 

.183 

196.19 

1.9816618 

1.681C619 

7.8116 

.217 

898.69 

1.9196618 

1.2896619 

6.9996 

.191 

167.19 

1.3316618 

1.1296619 

6.1796 
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Table  le.  Atmospheric  Profiles  of  Temperature  and  Composition 
Derived  from  Rawlnsonde  and  Rocketsonde  Data 


EL 

PRSO,  TX 

PRESSURE 

TENPEROTURE 

HOLECULES/SO 

C7 

tH9l 

IK) 

C02 

N20 

03 

3SI.000 

307.05 

2.O90E«21 

2.079E«ia 

5.130E417 

311. OtO 

339.09 

2.100E«21 

1.7a2E«10 

5.767Eei7 

Esi.cai 

229.39 

1.750E*21 

1.0a9E«^18 

6.093Eel7 

Eoa.ota 

319.69 

l.O00E»21 

i.iaaE«ia 

7.20OE417 

iva.oaa 

311.09 

1.190E*21 

i.oioE^ia 

7,728E417 

laa.oaa 

303.99 

l.ei5E»21 

a.613E»17 

(.313E417 

isa.aaa 

300.05 

9.100E*20 

7.722€*17 

6.771E*17 

lao.oaa 

(90.05 

7.000E»20 

5.900€*17 

1.033E416 

ss.aao 

196.69 

6.690E«20 

9.603E»17 

1.070E*18 

aa.ooa 

(90.09 

6.300E«20 

5.306E^17 

1.119E«ia 

as.oaa 

300.69 

9.990E»20 

9.009E»17 

l.l77Eel8 

ea.aao 

303.05 

9.600Ee20 

0.752Eei7 

1.291Eel8 

fs.aae 

205.15 

9.290E«20 

0.055E»17 

1.392E«18 

Ec.oaa 

307.39 

O.C00E«20 

0.15aE^17 

1.071E«18 

as.aaa 

309.29 

O.590E«20 

3.a61E>17 

!.626E^18 

aa.oaa 

311.35 

0.200E720 

3.560Eei7 

1.817Eel8 

ss.oaa 

313.65 

3.e9OE«30 

3.267E»17 

2.032E418 

ac.oaa 

315.95 

3.500E«20 

3.970E*17 

2.287E*ia 

ks.ooa 

316.29 

3.190E420 

3.673E»17 

2.979E418 

aa.oao 

316.09 

2.600E»20 

3.376E«17 

2.098Eei8 

36.aaa 

317.35 

2.530E»20 

2.13aEei7 

3.162Eel8 

33.000 

317.99 

2.310E«20 

1.960E«17 

3.011E«18 

30.000 

310.09 

2.100£«20 

1.7a2£^17 

3.668£*18 

sa.ooo 

330.29 

l.t20E*<20 

1.500E»17 

0.867E4ia 

23.000 

222.15 

1.610E*20 

1.366E»17 

0.012Eel8 

20.000 

223.65 

1.000E»30 

i.iaaE»i7 

0.707E418 

10.000 

325.65 

1.260E«20 

1.0E9£«17 

!.092E«18 

IE. 000 

320.19 

1.120E«20 

9.90034^16 

5.320E«ia 

10.300 

330.19 

l.C01E«20 

8.090E«16 

!,577E*18 

12.000 

330.39 

6.000Eei9 

7.120Eei6 

5.931E418 

10.000 

230.39 

7.COOE«19 

9.900E416 

C.206£«18 

a.ioT 

337.95 

5,703E*19 

0.039E416 

6.S90E«ia 

6.197 

203.35 

0.310E*19 

3.657E«16 

(.901E«18 

0.091 

350.09 

2.660E419 

2.O30Eei6 

2.398Eel8 

2.137 

260.69 

1.096E»19 

1.269E416 

1.138E417 

1.006 

375.35 

7.C02E«10 

9.976E«19 

3.071E4l( 

.930 

263.09 

3.76'SE«ie 

3.196E419 

1.079E416 

.322 

390.09 

2.250E«lfl 

1.913E419 

l.200Eel5 

.203 

296.19 

l.SOlEeia 

1.6I1E419 

7.011E419 

.217 

290.09 

1.519E*16 

1.2a9E«19 

0.999EelS 

.190 

207,15 

1.330Eeie 

1.129E419 

*.175E415 
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Table  It.  Atmospheric  Profiles  of  Temperature  and  Composltloa 
Derived  from  Rawlnsonde  and  Rocketaonde  Data 


AlOUQUCOQUCf  NN 

PRCSStME 

TCOPCOOTMC 

NLECULCS/$0 

CO 

fN8) 

<K} 

C02 

N20 

03 

1 

3ff.MI 

269.09 

2.690Ee21 

2.079EelO 

9.130Eel7 

911. IN 

130.69 

l.lOOCail 

1.702Eel0 

!.767Eel7 

ZN.IN 

029.99 

1.790E»21 

1.609Eei0 

C.69SE*17 

IN. IN 

021.19 

1.600E421 

l.lOOEelO 

7.206E«17 

179. IN 

216.39 

1.229Ee21 

l.OOOEelO 

7.693Eel7 

191. ON 

207.09 

1.090E«21 

0.910Eai7 

f.l00E«17 

1II.N8 

200.29 

e.960Ee20 

7.267Eei7 

9.073Eel7 

IN. NO 

200.19 

7.000Ea20 

9.960Eei7 

1.093E«10 

91.900 

190.19 

6.679Ce20 

9.699Eai7 

1.097Eel0 

10.900 

202.09 

6.199Ee20 

9.297E^17 

1.136EeiO 

03.900 

100.09 

9.069Ee20 

6.960C»17 

l.lOOEeiO 

CO. 000 

201.19 

9.600Ee20 

6.792EV17 

1.29icn0 

79.000 

206.39 

9.290Ea20 

6.699Eel7 

l.S92EnO 

70.000 

207.29 

6.900E«i20 

6.190Eal7 

1.671EelO 

C9.000 

210.29 

6.390E«20 

3.061E^17 

1.626EaiO 

Cl. 900 

213.19 

6.309E»20 

3.693Eai7 

i.rooEvio 

CO. 000 

219.29 

6.200E«20 

3.966Eai7 

1.017E«10 

99.000 

213.99 

3.C90EN0 

S.267Eel7 

{.032EnO 

90.000 

216.29 

3.90IE«*20 

2.970E«17 

{.207Eai0 

%9.000 

219.99 

3.t90E»20 

2.673Eei7 

2.979EM0 

«0.000 

210.39 

2.C00Ea20 

2.376Eai7 

2.090E«10 

90.900 

210.99 

2.999E*20 

2.1C0Eai7 

3.169E«10 

31.300 

217.19 

2.'261E^20 

1.919E«17 

9.669EM0 

90.000 

220.69 

2.100E*20 

1.702Eai7 

9.E60Eei9 

26.000 

222.29 

1.020E^20 

1.966Eel7 

6.067Eel0 

23.000 

223.69 

l.ClOEaCO 

1.366E«17 

6.612E«10 

20.000 

229.39 

1.600E420 

1.10OE«i7 

6.707E«10 

10.000 

226.99 

1.260Ce20 

1.069E»17 

9.092EaiO 

10.000 

220.09 

1.120E*20 

9.906Eei6 

9.326Ene 

16.300 

230.09 

l.OOlEeOO 

0.696E«16 

!.977Eaie 

12.000 

230.09 

0.600Eal9 

7.120E^16 

9.931E010 

10.000 

230.09 

7.000E«19 

9.960E«16 

C.266E«16 

0.167 

237.99 

9.703Eel9 

6.039Eai6 

C.996EnO 

6.197 

262.39 

6.910Eel9 

3.697Eai6 

C.901E410 

6.091 

290.09 

2.066Eel9 

2.630Eai6 

9.390E«10 

2.137 

266.69 

1.696E«19 

1.269Eeie 

1.120EN7 

1.006 

279.39 

7.062Cei0 

9.976eei9 

3.071E«1E 

.930 

263.69 

3.766EelO 

3.196Eel9 

1.679Eel6 

.321 

290.09 

2.296E«ie 

1.913Eat9 

0.200En9 

.203 

296.19 

l.OOlEelO 

l.OOlEalS 

7.011Eal9 

.217 

290.69 

1.919E«10 

1.209Eal9 

6. 9992019 

.190 

267.19 

1.330Eei0 

1.129£«19 

6.179E019 

where 


is  the  radiance  received  at  the  sensor; 

V  is  the  mean  frequency  of  the  filter; 

8(7,  T)  is  the  Planck  Function  associated  with  frequency  7  and  temperature  T: 

▼(T)  is  the  transmittance  from  points  in  the  atmosphere  (where  the  temperature 
is  T)  to  the  point  of  observation.  This  transmittance  is  a  wei;;hted  aver¬ 
age  over  the  filter  function  shown  in  Figure  3. 

Initially,  the  four  rawinsonde  data  sets  obtained  at  White  Sands.  New  Mexico 
were  used  directly  in  a  numerical  analogue  of  Eq.  (1).  The  transmission  results 
were  computed  with  the  AFGL  Atmospheric  Absorption  Line  Parameters  Com¬ 
pilation  and  a  line-by-line  calculation  scheme  similar  to  the  one  described  by 
McClatchey  et  al^  but  applied  to  a  multi-layered  atmosphere  was  used.  The  trans¬ 
mission  results  were  calculated  in  this  way  from  each  pressure  level  specified  in 
the  rawinsonde  data  of  Table  1  to  each  lower  pressure  level  in  the  saime  sounding 
along  a  20*^  elevation  angle. 

It  was  felt  that  the  irregular  spacing  of  the  rawinsonde  pressure  levels  might 
lead  to  computational  errors.  This  could  be  a  particular  problem  if  large  inter¬ 
vals  occur  at  altitudes  where  the  transmission  changes  rapidly  from  near  unity  to 
near  zero.  In  order  to  avoid  this  problem  and  reduce  the  risk  of  a  substantial 
computational  uncertainty,  it  was  decided  to  define  a  much  finer  pressure  mesh 
and  then  to  interpolate  the  transmission  and  temperature  values  associated  with 
an  individual  rawinsonde  run  to  this  revised  pressure  mesh. 

Under  the  assumption  that  processes  in  the  atmosphere  are  adiabatic  (a 
reasonable  assumption  in  the  bulk  of  the  troposphere),  the  relation  expressed  as 

e 

Eq.  (2)  can  be  derived  (see  Gordon  )  where  k  is  the  ratio  of  the  universal  gas  con¬ 
stant  to  the  specific  heat  of  dry  air  at  constant  pressure. 

T  =  Const.  X  P"  (2) 

The  value  of  k  is  0.286  «  2/7.  To  the  extent  that  this  relation  is  valid,  it  can  be 
seen  that  equal  increments  of  p^^'^  would  yield  equal  temperature  increments.  An 
examination  of  the  hydrostatic  equation  indicates  that  equal  height  increments  will 
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also  result.  As  temperature  is  one  of  the  important  atmospheric  parameters 
governing  the  emitted  radiance,  we  have  chosen  equal  increments  of  p^  ^  as  the 
basis  of  increasing  the  number  of  atmospheric  levels  for  computational  purposes. 
In  this  way  we  have  increased  the  number  of  levels  from  a  typical  rawinsonde  set 
of  about  40  to  100  levels. 

Having  extended  the  independent  variable  array  as  indicated  above,  the  follow¬ 
ing  interpolation  procedures  were  devised  for  the  transmission  and  temperature 
respectively;  1)  The  transmission  computed  at  the  initial  rawinsonde  levels  was 
interpolated  exponentially  with  pressure;  and  2)  The  temperature  interpolation 
was  carried  out  on  the  basis  of  Eq,  (2),  but  with  the  exponent  determined  from  the 
pressure  and  temperature  values  in  adjacent  rawinsonde  levels.  Test  results 
were  computed  with  the  atmosphere  above  the  350  mb  level  divided  into  100  and 
200  levels  and  the  results  compared  with  those  based  on  the  original  rawinsonde 
set  of  about  40.  Changes  of  a  few  percent  occurred  when  the  100-level  results 
were  compared  with  the  40-level  results.  The  further  increase  to  200-levels 
caused  a  negligible  further  change  in  the  resulting  radiances.  Therefore,  all 
remaining  calculations  were  performed  with  100-level  atmospheric  models. 

Figures  9a  and  9b  provide  a  composite  of  the  calculated  radiances  based  on 
the  four  White  Sands  data  sets  and  the  direct  comparison  of  calculations  with 
measurements.  W'e  have  plotted  the  radiances  based  on  the  "local  temperature" 
(assuming  that  the  Planck  function,  B(v,T)  associated  with  the  local  temperature 
represents  the  total  radiance).  This  will  be  the  case  if  the  atmosphere  becomes 
opaque  in  a  depth  over  which  the  temperature  does  not  change  from  the  local 
value.  We  have  also  applied  Eq.  (1)  to  the  calculation  of  radiance  based  on  the 
rawinsonde/rocketsonde  profiles  obtained  at  each  of  the  indicated  times.  Although 
the  measured  radiances  are  seen  to  agree  rather  well  with  the  "local  temperature" 
radiance  values  in  the  range  from  290-150  mb,  it  can  be  seen  that  there  is  a  clear 
systematic  discrepancy  between  the  measurements  and  the  (presumably  more 
accurate)  radiative  transfer  calculations  (symbols  in  Figures  9a  and  9b),  The 
"local  temperature”  results  clearly  fail  for  pressures  less  than  125  mb,  but  we 
would  assume  our  radiative  transfer  calculations  to  represent  the  measured  radi¬ 
ances  rather  well.  The  disagreement  between  measured  and  calculated  radiances 
is  seen  to  increase  to  a  maximum  value  near  50  mb  of  about  20  percent.  It  is 
interesting  to  note  that  the  general  shape  of  the  calculated  and  measured  radiances 
match  even  though  their  absolute  values  do  not.  It  should  also  be  noted  that  the 
spread  in  calculated  radiance  values  associated  with  these  four  different  atmos¬ 
pheric  profiles  (see  Table  la  to  10,  is  small.  Thus,  it  seems  unlikely  that  some 
local  temperature  anomaly  associated  with  the  radiometer  flight  path  is  the  cause 
of  this  discrepancy. 


RADIANCE  PROFILE  (15/Afn) 
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Figure  9.  Computed  vs  Measured  Radiance  Profiles.  Computed 
points  are  the  results  of  application  of  the  radiative  transfer  equa¬ 
tion.  Computed  curves  result  from  the  equivalent  radiances  de¬ 
rived  from  local  temperature  values.  Results  are  provided  for 
the  four  different  White  Sands  Missile  Range  rawinsonde  data  sets 
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RADIANCE  PROFILE  (15/im) 
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Figure  10.  Computed  vs  Measured  Radiance  Profiles. 
Computed  points  are  the  result  of  application  of  the  rad¬ 
iative  transfer  equation.  Computed  curves  result  from 
the  equivalent  radiances  derived  from  local  temperature 
values.  Results  are  provided  for  the  Albuquerque  and 
El  Paso  raw  insonde  data  sets 


In  order  to  further  investigate  the  possibility  of  this  discrepancy  being  re¬ 
lated  to  an  atmospheric  non-uniform ity,  rawinsonde  data  from  two  nearby  stations 
(Albuquerque  and  El  Paso)  were  used  as  the  basis  of  radiance  calculations  again 
using  Eq.  (1).  The  results  can  be  seen  in  Figures  10a  and  10b  to  be  similar  to 
those  obtained  with  the  rawinsonde  data  from  White  Sands  Missile  Range. 


4.  CONCLUSIONS 

3 

In  view  of  the  previous  study  carried  out  by  McClatchey,  an  effort  has  been 
made  here  to  investigate  a  possible  source  of  the  identified  discrepancy  between 
computed  and  satellite-measured  radiances.  We  seem  to  have  generated  more 
questions  than  we  have  answered.  Although  the  calculations  in  the  range  from 
290-150  mb  only  differ  from  the  measurements  by  a  few  percent,  the  fact  that  the 
discrepancy  is  systematic  and  similar  for  all  rawinsonde  data  gives  cause  for 
concern.  At  higher  altitudes,  this  discrepancy  increases  to  intolerable  levels. 
How  can  we  hope  to  adequately  solve  the  inverse  problem  of  determining  the 
temperature  profile  from  radiometric  measurements  when  we  cannot  even  match 
our  calculations  and  measurements  under  conditions  of  a  well  documented  atmos¬ 
phere,  well-calibrated  instrumentation,  and  carefully  constructed  computation 
techniques  ? 

Let  us  look  at  the  implication  of  this  discrepancy  as  though  it  were  attributed 
to  transmission  errors  alone.  For  this  purpose,  we  might  inspect  the  results  in 
the  two  different  regimes  indicated  on  Figures  9  and  10.  Let  us  consider  sepa¬ 
rately  the  250  mb  and  50  mb  results: 

i)  At  250  mb 

In  view  of  the  general  agreement  between  the  "local  temperature" 
radiances  and  the  measured  radiances  in  the  region  from  290  mb- 150  mb,  it  is 
clear  that  an  extremely  large  transmission  error  would  be  required  to  explain 
the  discrepancy  between  measured  radiances  and  those  calculated  based  on 
Eq.  (1).  That  is,  the  atmospheric  distance  over  which  the  transmission  falls 
from  1.0  to  near  zero  would  have  to  be  represented  by  a  uniform  temperature. 
The  lapse  rate  as  indicated  by  the  rawinsonde  profiles  (see  Table  1)  is  rather 
steep  in  this  region.  An  error  of  50  to  100  percent  in  calculated  transmission 
would  be  required  to  drive  the  computed  radiances  into  agreement  with  the 
measurements. 

ii)  At  50  mb 

2 

The  measured  radiance  is  about  167  jiw/cm  sr  pm  compared  with  a 
2 

calculated  value  of  137.4pw/cm  srpm.  Our  calculations  clearly  show  that  the 
transmission  to  space  along  a  path  at  a  20°  elevation  angle  should  be  0.  31.  In 
other  words,  at  this  altitude  and  for  the  path  of  observation,  the  atmosphere  is 
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clt'nrly  not  opaque.  Let  uh  e.sfliiiate  the  net  trananil.'4Sit)n  inodincation  required 
ftir  the  raU'ulation  to  agree  with  the  inea.sureinent .  To  do  this,  we  assume  that 
we  ran  define  a  mean  Planck  radiance  for  the  entire  path  fi-om  sensor  to  space. 

We  would  have  (from  Kq.  (O),  Kr)  -  H(r,  T)(l  -  with  1(7)  137.4  and 

t(7)  0.31.  Thus,  the  mean  Planck  radiance  is  199.  1.  The  radiance  discrepancy 

would  imply  that  the  transmission  satisfy  the  equation,  197  199.  13  (1  -  t(7)). 

The  resulting  transmission  would  be  0.  Ui,  a  factor  of  2  lower  than  the  calculated 
value.  A  careful  analysis  of  our  transmission  calculation  accuracy  is  being 
conducted,  and  at  the  moment,  we  feel  that  this  kind  of  error  is  far  outside  the 
range  of  possible  transmission  uncertainty. 

In  both  the  250  mb  and  the  50  mb  measurements,  the  results  indicate  that 
the  atmosptiere  is  apparently  more  opaque  than  we  calculate.  An  alternative 
description  is  that  our  calculation  underestimates  the  near  field  radiances. 

These  results  could  also  be  explained  by  atmospheric  temperatures  which  are,  in 
fact,  higher  than  we  think.  It  could  also  result  from  some  additional  source  of 
atmospheric  radiation. 

It  IS  useful  to  compare  the  results  obtained  from  these  in-situ  measurements 
witfi  the  earlier  study  of  satellite-observed  radiances.  At  first  glance,  the  re¬ 
mits  might  appear  to  opfKise  <*ach  other,  the  satellite  study  yielding  a  positive 
calculation-measurement  discrepancy  and  the  current  study  yielding  a  negative 
calculation-measurement  result.  However,  a  common  thread  can  be  found.  The 
satellite  discrepancy  occurs  primarily  in  the  wing  of  the  carbon  dioxide  band 
where  the  satellite  sensor  channels  sound  the  troposphere.  An  effect  which  mani¬ 
fests  itself  as  an  apparent  increase  in  atmospheric  opacity  over  our  calculations 
would  tend  toward  observing  the  atmosphere  at  higher  altitude  where  temperatures 
(and  emitted  radiances)  are  less.  Both  these  discrepant  results  are  compatible 
with  the  interpretation  that  the  calculations  underestimate  the  near  field  contribu¬ 
ting  to  the  radiance  and  conversely  overestimate  the  radiance  from  the  far  field. 

In  summary,  we  feel  that  the  radiance  discrepancies  described  in  this  report 
are  real.  They  are  consistent  with  the  results  described  earlier  (Hef.  3)  In  con¬ 
nection  with  the  DM.SP  satellite  radiances.  P'urthermore,  we  feel  that  the  dis¬ 
crepancies  cannot  be  attributed  to  the  more  obvious  effects  of  Instrument  calibra¬ 
tion,  atmospheric  transmission,  or  lack  of  atmospheric  temperature  uniformity. 
An  examination  of  Kq.  (1)  leads  us  to  the  conclusion  that  the  only  remaining  source 
of  error  may  be  the  assumption  that  the  source  function  la  equal  to  the  Planck 
intensity.  This  possibility  has  sufficiently  large  ramifications  to  remote  sounding 
that  we  feel  it  should  be  investigated  further. 

We  intend  to  pursue  this  problem  further  with  the  launch  of  two  radiometers 
(one  operating  in  the  15  pm  and  one  In  the  4.  3  pm  COj  bands)  as  part  of  a  balloon 
flight  payload.  We  are  hopehil  to  make  measurements  starting  at  lower  altitudes 
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and  to  reduce  the  elevation  angle  of  those  measurements.  Furthermore  the  In¬ 
creased  opacity  of  the  4.3  jim  COg  band  Increases  the  altitude  at  which  our 
radiometer  can  be  used  to  measure  the  local  atmospheric  temperature. 
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